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We will summarize the progress in understanding the changes in the vector meson 
spectral density in nuclear medium using the constraint equations obtained from the Borel 
transformed dispersion relation and QCD Operator Product Expansion. We will discuss 
O ■ the results for the scalar mass shift and dispersion effects (three momentum dependence) 



for the light quark system {p,uj), the strange quark system (0) and the heavy quark 



O ! system (J/ip) in nuclear medium. For the light quark systems, a nontrivial change in the 
mass and width are expected, while the dispersion effects are found to be small. Existing 
model calculations for the dispersion effects are compared to the constraint equation in 
detail. Very small, but accurate mass shift is obtained for the heavy quark system. 



Q . 

^ . 1. Introduction 

The properties of vector meson in nuclear medium have been the focus of current in- 
^ ! terest due to their potential role to provide one with a direct observable of the nuclear 
^ I medium effects, associated with chiral symmetry restoration, via dileptons in p-A or A-A 
reactions Indeed dileptons from Relativistic Heavy Ion Collisions (RHIC)p[ seemed 
to suggest a non-trivial change of the vector meson spectral density in a hot/dense envi- 
ronment, which can be understood in terms of model calculations based on decreasing 
vector meson masses in hot/dense medium[|,^. However, some model calculations 
based on changes of the vector meson spectral densities obtained using effective hadronic 
model calculations with no decrease in the mass also seem to explain the main features 
of the CERES data. In all of the approaches, the central question is, how the spectral 
density changes in hot /dense matter P]. In this talk, I will discuss the results for the scalar 
mass shift and dispersion effects (three momentum dependence) of the peak position of 
the spectral density for the light quark system (p, u) , the strange quark system (0) and 
the heavy quark system (J/ip) in nuclear medium. 

1.1. Vector mesons in vacuum 

The distinctive features appearing in the dilepton spectrum in p-A or A-A reactions 
are the vector meson peaks; the p, (p and the J/ip, whose (mass, width) in the vacuum 
are (770, 150), (1020, 4.4), (3100, 0.086) MeV. It is interesting to compare the phase 
spaces of their decay into their corresponding two pseudo particles. For the p, its two 
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pion decay, which accounts for most of the total width, has a large phase space. This is 
so because the pion is a Goldstone boson and has a small mass. For the J/ip, its decay 
into D-mesons are forbidden, because the D mesons are not Goldstone bosons and two 
times its mass is greater than the mass of J/ip. For the 0, the situation is something in 
between and its decay into two kaons has very little phase space. Hence, chiral symmetry 
breaking is partly responsible for the large difference in their width. As for the masses 
of the vector mesons, among other models, QCD sum rules provide an indirect relation 
to QCD condensates. For the light quark system, (gg) is dominantly responsible for its 
mass, for the strange quark system, {ss) is responsible, and for the J/V', the charmed 
quark mass and a small non perturbative contribution from the gluon condensate (— G^) 
are responsible. 

Therefore, if chiral symmetry gets restored at finite density or temperature, non-trivial 
changes will occur to the masses and widths of the vector mesons. 

1.2. Quark condensate at finite density 

The temperature dependence of the quark condensates have been calculated long ago 
in the lattice |l9|,[T0[|. The result is that the light quark condensate will go to zero above the 



critical temperature. For heavier quark masses, the changes become smaller. At present, 
due to technical reasons, no lattice result exist at finite density. In all what follows, we 
will use linear density approximation 

(O)p„ = (0)o + p„x(0);v (1) 

where O is any operator, p„ the nucleon density, and the subscripts pn,0,N denotes the 
nuclear, the vacuum and the nucleon expectation values. Then, we have the following 
model independent result |TllJT^ for the quark condensate 



{<iq)p. = {m)o + ^Pn ~ (gg)o ■ (i - 0.2 (2) 
and for the strange part, 

{ss)p. = {ss)o + y^Pn ~ {-ss)o ■ (l - 0.05 ^) , (3) 

where, y = 2{ss)n/{{uu)n- + {dd)^) and po is the nuclear matter saturation density. One 
notes that already at nuclear matter, we have partial restoration of chiral symmetry; 
namely, the chiral order parameter is reduced by 20%. Nuclear matter provides a stable 
environment with non trivial vacuum changes. Hence, if anything happens to the vector 
mesons at high temperature or density, the tendencies will already be apparent at nuclear 
matter. 

1.3. Gluon condensate at finite density 

The temperature dependence of the gluon condensates has also been calculated on the 
lattice |ll3HT6|. The result is that it will stay almost constant up to the critical temperature 



and then reduce to about 60% of its vacuum value above the critical temperature. The 
leading density behavior can be obtained from the trace anomaly relation to leading order 
in a. 



Tli = -l^G' + j:mm (4) 
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Taking the nucleon expectation value of it and using the most recent determination of 
the nucleon mass in the chiral limit ITTIl, we have 



(-GV = Ag^-(i-0.05^) (5) 

Hence, we have a non-trivial change in the gluon condensate, although its relative change 
is smaller than the case of the quark condensate. 

As we have seen, the condensates have non-trivial change in the nuclear medium, which 
will have a non-trivial effect on the vector meson properties in nuclear medium. Exper- 
imentally, there are attempts to produce and observe the decay of vector mesons inside 
a heavy nuclei and also produce meson bound nuclei. Hence, the theoretical study and 
the subsequent experimental verification of changes of vector meson properties in nuclear 
medium will provide a solid basis for understanding the hadronic effects in RHIC and 
medium effects in general. 

2. QCD constraints 

There have been many model calculations to study vector meson properties in nuclear 
medium. Here, we will avoid any model calculation and derive a constraint equation that 
any model calculation should satisfy. The foundation of this approach was laid in ref. 

Consider the correlation function of the vector current = g7^g at finite density; 

H,.(a;2,q2) = ijd'xe^%T[J,{x)JM])p^- (6) 

Here q = {uj,q). In what follows, when we give result for explicit vector meson, we will 
use the currents JP'"^ = ^{u'j^u =F dlfid) for the p,u mesons, = S7^s for the and 
J^/'^ = C7^c for the J/ip. 
In general, because the vector current is conserved, the polarization tensor in eq.(BD 



will have only two invariant functions jTS 



U^^iu', q2) = H^g^Pj, + H^g^Pj,, (7) 

where we assume the ground state to be at rest, such that, Pqq = Pq^ = P^ = 0, 
^5 = ~ ^i^j/q^ and P^^ = {q^,qu/q^ - Of^u - P^^)- When q ^ 0, H^ = Hy, as in the 
vacuum. 

We will make a small q expansion of the correlation function and look at its energy 
dispersion relation at fixed q, 

ReHL,T(c^^ q^) = Re (n°(cj2, 0) + H^ ^.(cj^ 0) q^ + 

T~2 ^ + q + ■■ ' (8) 

where p('U^, q^) = l/7rImH^(M^, q^), and R denotes the retarded correlation function. We 
will construct a constraint equation for H'',n];, and H;^. For H]^,Hy, we will only look 



at the "non-trivial" q dependence ||19[|. A simple method to extract the "non-trivial" q 
dependence is to express the polarization function in terms of Q^,q^ (H(Q^,q^)) and 
extract the linear q^ term. 
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In general for eacli polarization functions (11°, H^), the OPE looks as follows, 

n(cu2) = 5]a(o„). (9) 

n 

Here the On are operators of (mass) dimension n, renormalized at a scale /i^, and Cn are 
the perturbative Wilson coefficients, which for the light quark system can be written as 
Cn = (_j2")n/2 and for heavy quark system as C„ = 

Let us first discuss the light quark system. After the Borel transformation, the disper- 
sion relation for any one of the polarization function(n°, U^p) becomes 

B.T. Ren(M2) = J dsp{s)e~''^'\ (10) 

The left hand side, which is the Borel transform (B.T.) of the OPE, looks like the following 
when including operators up to dimension 6, 

B.T. Ren(M2) = B.T. Ii{M\^rt + + ^(Oe) (11) 

The truncation is valid as long as is sufficiently large. The minimum M^^^ is usually 
determined by requiring the correction from higher dimensional operators to be less than 
30% of the perturbative contribution. Now the constraints for the spectral density would 
be eq.(|TO|), applied above > M^^^. 

As can be seen in eq.(|lO[), the Borel transformation also changes the weighting factor 
of the spectral density to an exp{—s/M'^). This has the following advantage for practical 
applications of our constraint. For small values of the Borel mass, the contribution of the 
spectral density at larger energy is exponentially suppressed. Consequently, in a model 
calculation, one can concentrate on the changes of the spectral density near the vector 
meson mass region and below and model the higher energy part with a simple pole like 
contribution. 

The constraint equation in the vacuum are well satisfied by the spectral density in the 
vacuum. As we will see, in most cases the changes in the operators (0„) are known. 
Hence, starting from the vacuum form of the spectral density (ImH), we can study what 
changes are consistent with the constraint equation. This provides model independent 
QCD constraints that any model calculation should satisfy. One can go one step further 
and try to parameterize the spectral density with a simple delta function type of pole and 
a continuum and determine the changes in the parameters. 

3. Light quark system p,uj 
3.1. Ho 

The operators that dominate the change in the OPE in eq.(0) for the light quark 
system (p, tu) are the quark operators]^ 

(O4) {qii^D^q) oc J dxx[q{x) + q{x)] 

(Oe) ^ {{qq?)<x{{qq)\ + 2pn{qq)o{qq)N (12) 



The first equation of eq.(|12D dominates the changes in the dimension 4 operators and is 
related to the well known second moment of the quark distribution function. The second 
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equation of eq.(|T2D, for which we have used the ground state saturation hypothesis[^, 
dominates the changes in the dimension 6 operators. 

Using a delta function assumption for the spectral density in the constraint equation 
gives the following result for the scalar mass shift at q = , 

""""^^"^ (0.16 ±0.06)^ (13) 



mvipn = 0) ■ ' Po 



This result is also consistent with other model calculations! 23- 25 1 or the Brown- Rho 
scaling argument [Q. 

A detailed comparison of the constraint equation in eq. ( p!OD with a hadronic calculation, 
based on chiral SU(3) dynamics with explicit vector mesons were performed in [^]. The 



result shows a very good agreement between the OPE and the phenomenological spectral 
density put into the constraint equation in eq.([10D. However, the hadronic calculation 
gave a large increase in width with a small decrease in mass for the p and a large decrease 
in mass with a small increase in width for the u. Hence only the result for the u is 
consistent with the result in eq.(0). Later it was found that this was a general result, 
given the uncertainty in the ground state hypothesis for the four quark condensate in the 
medium[^; namely, that there exists a band in the mass vs. width plane that satisfies 
the constraint equation in eq. (pUj) . 

3.2. Ht 

The constraint for the nontrivial q dependence in eg. (p!0D has no H(M^)pert and has con- 
tributions from operators with explicit spin index [|r^. The contributions from dimension 
4 operators are related to the twist-2 matrix elements and are well known. The dimen- 
sion 6 operators are dominated also by the twist-2 matrix elements. Hence the constraint 
equation has little uncertainty coming from the OPE. The operators that dominate are 



(O4) {qi^D^q)N oc J dxx[q{x) + q{x)] 

(Oe) ^ {q^^D^D^D0)N (X J dxx^[q{x) + q{x)] (14) 

Using a delta function assumption for the spectral density and allowing the param- 
eters to change to leading order in density and in q^, we find the following non-trivial 
momentum dependence in the peak position! 



""^^^"^ , = 1 - (0.023 ± 0.007) f-S-V ^ 
mp{pn = 0) V0.5/ Po 



(0.016 ±0.005) (^) — , (15) 



m^(Pn = 0) ^ V — ^ V0.5y Po 

where q is in GeV/c unit. This shows a very small momentum dependence compared to 
the expected scalar mass shift in eq.(0). 

A detailed comparison of the constraint equation in eq. ([T0|) for the momentum depen- 
dence for the transverse direction has been made [ 27 1 to the hadronic calculation, where 



the vector-meson nucleon scattering amplitude is obtained by resonance saturation in the 
s-channel. The result shows that the existing model calculations tend to overestimate the 
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constraint. This is due to the large p — N — A(1232) couphng, which is obtained from 
the Bonn potential |^8[. However, the existing calculations used a non-covariant monopole 
form factor |^,^| normalized off shell 



F(q) 



A2 



q 



2 ■ 



(16) 



On the other hand, the large p — N — A(1232) coupling in the Bonn potential is defined 
with a dipole form factor normalized at the on shell point of the vector meson. 



2\ 2 



(17) 



This reduces the Delta contribution to the rho meson self energy at the invariant mass 
around — m-^v by approximately a factor of 4. After this correction, we find that the 



model calculations give very good agreement with the constraint equation |27] 



3.3. Ul 

The constraint for the longitudinal direction is dominated by twist-2 quark and gluon 
operators {q-f^,D^^..D^^q)N oc /c/xx"-^[g(x) + q{x)], {G'^^Df,^..G^^a)N oc J (ixx'^-^gix), 

Using a delta function assumption for the spectral density and allowing the parameters 
to change to leading order in density and in q^, we find||19||. 



, = 1 - (0.004 ± 0.002) f A V ^, 



:i8) 



for both the p and u. This is a very small effect and no detailed comparison with any 
hadronic calculations exits yet. 

4. Strange quark system 

The OPE in the constraint equation for the meson is dominated by (rrisss) for the 
scalar mass shift and {s'-y^Di.s) for the momentum dependence. Assuming a delta function 
ansatz for the pole, we find 



""'^^^"^ , = 1 - (0.03 ± 0.015)^ 
m4,[pn = 0) Po 



(0.0005 ±0.0002) 



q 



Pn 



0.57 Po 



(19) 



for the transverse vector meson and the momentum dependence for the longitudinal is 
about a factor of two larger. 

5. Heavy quark system J/-ip 



For the heavy quark system, we look at the constraints from the moments p0|j31 



d 



Pis) 



ds. 



(20) 
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The changes in OPE in nuclear medium for M„ is dominated by the change in the gluon 
condensate in eq.(^). 

To study the J / ip at rest in the nuclear matter, we will approximate the spectral density 
with a delta function for the lowest state. This is valid even in nuclear matter, because 
for a J/ip at rest, inelastic interactions with nucleons such as J/ip + N Z) + Ac do 
not occur. With this assumption, the constraint equation allows for the determination of 



mass shift of the J/^p in nuclear matter. We find|^ for the mass shift 

Amj/^~-7MeV. (21) 

This corresponds to small J/ip- and r^c-nucleon scattering lengths a = — /irAm/(27rpiv) — 
(0.1 — 0.2) fm [fir is the meson-nucleon reduced mass). Our results for the mass shifts of 
the lowest cc states are surprisingly close to those reported in ref. |l33|-p6|. 



Although the expected mass shift is small, the result has little uncertainty coming from 
OPE and puts reliable constraint on charmonium mass shift which should be met by 
further studies of heavy quark systems in dense matter. 



6. Conclusion 

We have derived and explored the consequences of model independent constraints for 
the vector meson polarization at q = and q 7^ for all the vector mesons p, u, 0, J/ip- 
Most of them have very little uncertainty in the OPE side of the constraint equation and 
can be used as reliable constraints on all model calculation of the vector meson properties 
in medium. 
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